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G. X. Peng 1,2,3 , X. J. Wen 2 , Y. D. Chen 2 

1 China Center of Advanced Science and Technology (World Lab.), Beijing 100080, China 

2 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100039, China 

3 Center for Theoretical Physics MIT, 77 Mass. Ave., Cambridge, MA 02139-4307, USA 

Recent publications rule out the negatively charged beta equilibrium strangelets in ordinary phase, 
and the color-flavor locked (CFL) strangelets are reported to be also positively charged. This letter 
presents new solutions to the system equations where CFL strangelets are slightly negatively charged. 
If the ratio of the square-root bag constant to the gap parameter is smaller than 170 MeV, the CFL 
strangelets are more stable than iron and the normal unpaired strangelets. For the same parameters, 
however, the positively charged CFL strangelets are more stable. 
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After the acceptance of quantum chromodynamics as 
the fundamental theory of strong interactions, it became 
extremely significant whether a deconfined phase of mat¬ 
ter consisting merely of quarks would be possible. The¬ 
oretical investigations show that strange quark matter 
(SQM), which is composed of u , d, and s quarks, might be 
absolutely stable HiS- Because small lumps of SQM, 
the so called strangelets, could be produced in modern 
relativistic heavy-ion collision experiments, their charge 
property has attracted a lot of interest Q. 

Originally, SQM is believed to show up with some small 
positive charge @|. In June 1997, however, Schaffner- 
Bielicli et a 1. demonstrated that strangelets are most 
likely heavily negatively charged In June 1999, it was 
shown that negative charge can lower the critical density 
of SQM J6j. In July 1999, Wilczek mentioned an “ice- 
9”-type transition |7], which was picked up by a British 
newspaper. Not long ago, in response to public concern, 
an expert committee published a report |8l, which got 
positive comments as well as criticisms 11 ( 1 . In fact, 
the strangelets in Ref. [|] are not in (3 equilibrium which 
drives the system to flavor equilibrium, and negatively 
charged strangelets in normal phase have been ruled out 
by a recent publication 112 - 

Much progress has been achieved recently by the in¬ 
troduction of color supperconductivity [lj, 113 ■ It has 
been shown that bulk SQM with color-flavor locking 
is electrically neutral El- Immediately, Madsen found 
a solution to the corresponding system equations of 
strangelets, where color-flavor locked strangelets are pos¬ 
itively charged 0, and they might be a candidate for 
cosmic rays beyond the GZK cutoff [Tg| . 

Very recently, it is shown that CFL phase can exist 
only when the ratio of the squared strange quark mass 
to chemical potential, i.e., m 2 //x, is smaller than a crit¬ 
ical value about 2 times the gap parameter A 0 . it 
is therefore of interest to study if there is some similar 
criterion for CFL strangelets given that surface effects 
become quite important. Fig. ^ explicitly shows the ra¬ 
tio for various parameters. The solid lines are for the 
CFL strangelets reported in Ref. 0. These strangelets 
are positively charged. At the same time, there are new 
solutions (the dashed lines) which are slightly negatively 
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FIG. 1: The ratio of squared strange quark mass to chemi¬ 
cal potential for different strangelets with various parameters. 
The solid lines are for the CFL strangelets in Ref. 0. The 
dashed lines give the new solutions of CFL strangelets re¬ 
ported in this letter. Unlike the previous strangelts, which 
are positively charged, these new strangelets are slightly neg¬ 
atively charged, or nearly neutral. Parameters are indicated 
as (A, B 1 / 4 ,m s ) in MeV. 


charged or nearly charge-neutral, to be discussed in de¬ 
tail below. It is found that the stability of the CFL 
strangelets can be judged by the ratio \f~B/ A, i.e., the 
square-root bag constant to the gap parameter. If the 
ratio is less than about 170 MeV, these CFL strangelets 
are more stable than 56 Fe, i.e., the energy per baryon is 
less than 930 MeV. 

As done in Ref. 0 , the thermodynamic potential den¬ 
sity is written as Ul = flf + fl pa ir + B. The paring con¬ 
tribution is flpair = —3A 2 /l 2 /7r 2 with A being the paring 
gap and (1 = (fx u + Pd + p s )/3 being the average chemical 
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potential of quarks. The normal quark contribution is 


ilf = 


J (y\Jp 2 + m i - Pi'j rl '{p- m i, R)dp (1) 


with the density of state n'(p , rrii, R) given in the multi¬ 
expansion approach 0 by 

m V 
p )\ ’ 
( 2 ) 

where g=6 is the degeneracy factor for quarks, and 
the functions /s 0 and fc 0] are given, respec¬ 
tively, by fs(x) = — arctan(a:)/(47r 2 ) and fc(x ) = [1 — 
(3/2)arctan(x)/x]/(127r 2 ). 

The common Fermi momentum v is a fictional inter¬ 
mediate parameter. It does not fully specify the quark 
number density, as it does in the unpaired case. The basic 
requirement is that it can not be negative. As a general 
practice, it is determined by minimizing f 1 at fixed radius 
R, i.e., 


n'(p, rrii, R) = g 


P 2 , 3 , (m 
2tt 2 + R fs V p 


P+^/o 


an 

dv 


n(v,m.i,R) 

i=u,d,s 


is 2 + mj - iii 


= 0. 


( 3 ) 


The number densities rii (i = u , d, s ) for quarks are 


rii = 


dfl 


dpi 


R 


dfl dfl dv 
dpi dv dpi 


( 4 ) 


Because of Eq. 0 , the second term vanishes, while the 
first term gives rii = gv 3 /(6n 2 ) + 3n i} s/R + &n- h c/R 2 + 
2A 2 p/tt 2 , where 


ni,S 


n i,c 


9mj 

8tt 2 
gm 
107T 2 


<t>i - tan (fn - - (j)^ tan 2 fa , 

4>i + ^ tan fa ~ <t>i) tan2 <t>i 


( 5 ) 

( 6 ) 


is the same as in Refs. 00 for bulk CFL quark matter. 
For a given baryon number A, one naturally has 



i=u,d,s 


3 A 
AnR 3 


( 8 ) 


To maintain mechanical equilibrium, one must require 
that the pressure is zero, i.e., 


— ?§- 


( 9 ) 


Please note, there is an extra term when it is compared 
to the normal case P = —12. This is because of the direct 
radius (or volume) dependence of the thermodynamic po¬ 
tential density. 



FIG. 2: Charge of strangelets. The horizontal axis is the 
baryon number A. The vertical axis is the electric charge Z to 
A 2 / 3 . There are three kinds of CFL strangelets marked with 
CFL slet-1 (dashed line), CFL slet-2 (solid line), and CFL 
slet-3 (dotted line). They are, respectively, charge-positive, 
negative, and nearly neutral. The ordinary strangelets are 
also plotted (dot-dashed line). 


with fa = arctan(i//mi). 

The chemical potentials pi and the radius R are the 
independent state variables. For a given baryon number 
A, one should give these quantity to fix a strangelet. To 
have chemical/weak equilibrium, the chemical potentials 
Pi satisfy pd = p s = Pu + Pe, maintained by reactions 
such asTt + d <-> s + w, u + e~ <-> d + v e . Because 
the strangelet radius is much smaller than the Comp- 
tom wave length of electrons, the electron’s number den¬ 
sity, and accordingly the chemical potential p e , must be 
zero. Therefore, strangelets in perfect (3 equilibrium al¬ 
ways have p = p u = pd = p a = p- Consequently, Eq. @ 
gives 

p = ^ ~2n'(v,mi,R ) ^Jv 2 + m 2 / '^2n\v,m i ,R). (7) 

i i 

When R —> oo and m u = m,d = 0, t his equatio n gives 
p = (2v+ \Jv 2 + m 2 )/3 or v = 2p— \Jp 2 + m 2 / 3, which 


For a given baryon number A , one can solve the three 
equations 0 , 0 , and 0 for p, v , and R. Then the 
overall electric charge is Z = V{2n u /3> — rid/?> — n s /3) 
with V = 4ttR 3 /3 being the volume. Numerical re¬ 
sults are given in Fig. 0 for parameters A = 150 MeV, 
B = (155 MeV) 4 , and m s = 150 MeV. It is found that 
there are three solutions for each given baryon number A. 
The strangelet corresponding to the first solution (dashed 
line) is positively charged. It is just the one that has been 
previously found in Ref. [TH . The strangelet correspond¬ 
ing to the second solution is negatively charged (solid 
line), and the third solution is nearly neutral (dotted 
line). For convenience, these three solutions are marked, 
respectively, with CFL slet-1, 2 and 3. The ordinary 
strangelets without color-flavor locking have also been 
plotted in the same figure for comparison purpose. It 
can be seen that the charge of the CFL slet-1 is approxi¬ 
mately proportional to A 2 / 3 , while that of the CFL slet-2 
or 3 is nearly proportional to A 1 / 3 . 
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To have a better understanding of the three solutions, 
let’s take some mathematical analysis. 

First, assume the common Fermi momentum v is much 
bigger than the strange quark mass, i.e., m s /v <C 1. In 
this case, one can take the limit of m Ut d — 1 ► 0 first, and 
then expand to a Taylor series with respect to m s on 
all the above expresions, to get simple expressions. The 
expansion of the pressure is 


gm s 

107T 2 

9" 


irg 4v(2 fj, — v) 

tf 2 
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87T 2 


R 


v 2 {4fi — 3z/) — 


3A 2 
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For the quark number densities, they are 

2A 2 gv | gv 3 
Uu ' d ~ “ ~4tt 2 R 2 + 6^ 


-B. (10) 
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and 
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( 12 ) 



0.37 


0 CFL slet-1 

0.36 

\ CFL slet-2 





0.35 

’ 

- 

0.34 



0.33 

- - 

- 


:/ <j 

^ (a) 

0.32 

(b) 

..1_l_i i mill..1.. 

0.31 

• .. . 


0 0.31 . 

10° 10 1 10 2 10 3 10 4 10° io 1 io 2 io 3 io 4 


A A 



It is obvious from Eqs. m and m that n s is smaller 
than n Ut d because n s has an extra negative term. The 
corresponding strangelet, CFL slet-1, is thus positively 
charged. 

Secondly, assume v is modest, i.e., it is smaller than 
m s but larger than m Ut d- In this case, one can still take 
the limit of m Ut d 0 for u/d quarks. But for s quarks, 
expressions should be expanded according to z/, rather 
than m s - Accordingly, Eq. © becomes 


P = 


3A 2 


gm s v 
07r 2 P 2 


(13) 


The u/d quark number density is still the same as Eq. 
m For s quarks, however, one now has 


FIG. 3: Quark fractions of different strangelets. Figures (a)- 
(c) are for the three kinds of CFL strangelets. Figure (d) is 
for the ordinary strangelets. The vertical axis for each figure 
is the quark number density in unit of the total quark number 
density, or the ratio of the corresponding quark number to the 
total quark number. 


Namely, the three flavors of quarks are nearly equal in 
this case. The corresponding strangelet, CFL slet-3, is 
almost neutral. 

Naturally, the above expanded expressions are merely 
approximate. Real calculations have been performed by 
directly solving the original system equations. Fig. [3] 
shows the quark fraction in different phases. They are 
qualitatively consistent with the above analysis. 


2A 2 3gv /8 ttv \ 
7T 2 ^ + 16nR 2 \3 m s ) 


3 gv 2 

107TP 


(14) 


Please note the curvature (the R~ 2 term) contribution. 
It is negative for u/d quarks in Eq. dill) . However, it is 
positive for s quarks in Eq. ra. This makes n a bigger 
than n U} d- Consequently, the corresponding strangelet, 
CFL slet-2, is negatively charged. 

Thirdly, assume v is extremely small so that expansion 
can be done with respect to v for all the three flavors. In 
this case, the pressure gives 

P = ^rV 2 + ^(»-m)^- B , (15) 

where fh = (m u + nid + m s )/ 3, while the quark number 
densities are 


Tl u ~ ~ n s 


2A 2 

—+ 

7 T 2 


gy 

27 t 2 R 2 ' 


( 10 ) 



FIG. 4: Parameters for CFL strangelets to be more stable 
than 56 Fe. (a) is for CFL slet-2 and 3 while (b) is for CFL 
slet-1. The full dot indicates the parameters in this paper. 





































4 


Now we discuss the determination of parameters. For 
the u/d quark mass, we take m u = 5 MeV and rrid = 10 
MeV, which are closer to the accepted current mass 
of light quarks ;22|. Decreasing the u/d quark mass 
has little effects on CFL slet-1 and 2, while the charge 
of CFL slet-3 becomes smaller and smaller until it is 
charge-neutral. The strange quark mass is expected to 
be density-dependent, lying between the current mass 
~ 100 MeV and the vacuum constituent quark mass ~ 
500 MeV. The A value varies from several tens to sev¬ 
eral hundreds of MeV in literature. For example, it can 
range from 20 MeV to 90 MeV {2^, or from 50 MeV 
to more than 100 MeV @]. Sophisticated treatments 
of the instanton interaction, including form factors from 
suitable Fourier transformation of instanton profiles, give 
larger values for A, as large as more than 200 MeV 
|25j. Therefore, we treat A as a free parameter in the 
present investigation. For the above calculations in Figs. 
2 and 3, we have taken A = 150 MeM i? 1 / 4 = 155 MeV 
(This B value was also used in Ref. [2l|), and m s = 150 
MeV. How these parameters influence the stability of 
CFL strangelets will be discussed a little later. 

Although the ‘common Fermi momentum’ v in CFL 
slet-2 and 3 is small, the chemical potential /i is still 
large. To get an approximate expression for /i from the 
equality P = 0, one can take v = 0 in Eq. o or G3. 
resulting 


7T \/~B 


(17) 


For the parameters chosen for Figs. 2 and 3, Eq. 03 
gives n ~ 290 MeV, very close to the actual value from 
the numerical calculation, which gives the ratio m 2 /fi to 
be about 77 MeV. On the other hand, /i varies in the 
range of 240 — 263 MeV for CFL slet-1. 

The radius of CFL slet-2 and 3 can be approximately 
expressed as 


-Rslet —2,3 


3V3A \ V3 

8A VB J 


(18) 


This equation means R oc A 1 / 3 , which is a known fact in 
nuclear physics, one may perhaps imagine A, composed 
of ( uds ), as the simplest CFL slet-3. The H particle |26i . 
composed of ( uuddss ), is probably the next simplest CFL 
slet-3. 

For information on the stability of CFL strangelets, 
we should investigate the energy per baryon E/n b- It 
is generally a function of A,m s , A, and B , i.e., E/n b = 
/(A, m s , A,B). If E/nb is less than 930 MeV (the mass 
of 56 Fe divided by 56), the strangelets are absolutely 
more stable than normal nuclear matter. Otherwise, they 
are meta-stable or unstable. The full line in Fig. Ha) 
gives A as a function of B at A = 20, m s = 150 MeV, 
E/nb = 930 MeV. In fact, this line does not depend 
strongly on the concrete values of m s and A. Because 
E = n + H + 3/zrib and « —P = 0, one has 


E/nb ~ 3/i. With a view to Eq. <11711 . we immediately 
have E/nb ~ V?>itVB/ A. Therefore, if 


A? < HA! * 170 MeV, 
A 7T 


(19) 


then the parameter pair (B, A) is located in the up-left 
part of Fig. Ha), and the new strangelets are more stable 
than iron. For CFL slet-1, a similar solid line is plotted 
in Fig. Hb). For different m s and A, this line moves a 
little up-left (bigger m s e.g. the dotted line for m s = 180 
MeV) or down-right (smaller m s ). However, the line in 
Fig. Ha) is always located in the region where CFL slet-1 
is more stable than 56 Fe for reasonable m s . Therefore, 
if the condition Eq. m is satisfied, all the three kinds 
of CFL strangelets are more stable than 56 Fe, and also 
more stable than the normal unpaired strangelets. As 
for the comparative stability between the three kinds of 
CFL strangelets, it depends on the pairing parameter A. 
If one uses the same A for all the three, then the slet-1 is 
more stable. In this case, however, the former is denser. 
Because investigations have shown that A depends on 
dendity, most probably increases with increasing densi¬ 
ties |25|, the comparative stability of the three kinds of 
CFL strangelets needs to be further studied in the future. 

CFL strangelets which are more stable than 56 Fe may 
have far-reaching consequencies. The slet-1 can provide 
an alternative explanation for cosmic rays beyond the 
GZK cutoff 16]. The slet-3 is nearly neutral, and so 
might be a candidate for the miracle dark matter in our 
universe. The slet-2 and 3 are more stable than the nor¬ 
mal unpaired strangelets, and so may have chances to be 
produced in the modern heavy ion collision experiments. 
However, they are unable to transform our planet into a 
strange star for the following two reasons. First, the pos¬ 
itively charged slet-1 is the energy minimum for the same 
parameters. And secondly, when the electron’s Compton 
wave length (ss 386 fm) is reached, the constraint n e = 0, 
(or, equivalently, = Hd = Us) is no longer valid, and 
so the strangelet will be neutralized and ceases to expand 
its size. 

It should be emphasized that the strangelets reported 
here are different from the previous ones fillfl5l in that 
their electric charge is opposite. The new strangelets 
are also different from the heavily negatively charged 
strangelets in Ref. @. There the strangelets were not in (3 
equilibrium and had no color-flavor locking. It was inves¬ 
tigated how the metastable candidates might look like if 
they are assumed to be stable against strong hadronic de¬ 
cay and subsequently against weak hadronic decay. Here 
the strangelets are assumed to be in perfect (3 equilibrium 
and considered as having the possibility of absolute sta¬ 
bility. However, the concrete values should not be taken 
seriously, and further studies are needed. 

In summary, there exist new solutions to the system 
equations where CFL strangelets are slightly negatively 
charged or nearly neutral. If the ratio of the squared bag 
constant to the gap parameter is smaller than 170 MeV, 
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CFL strangelets are more stable than the normal nuclear 
matter and ordinary unpaired strangelets. 
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